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C/3 , Abstract. Red supergiants (RSGs) are a He-burning phase in the evolution 

of moderately high mass stars (10-25 solar masses). The evolution of these 
stars, particularly at low metallicities, is still poorly understood. The latest- 
, type RSGs in the Magellanic Clouds are cooler than the current evolutionary 

J> ■ tracks allow, occupying the region to the right of the Hayashi limit where stars 

■^j- ' are no longer in hydrostatic equilibrium. We have discovered four Magellanic 

t-H , Cloud RSGs in this region that display remarkably similar unusual behavior. 

0^ ■ All of the, show considerable variations in their V magnitudes and effective 

CN | temperatures (and spectral types). Two of these stars, HV 11423 and [M2002] 

qq , SMC 055188, have been observed in an M4.5 I state, considerably later and cooler 

■ than any other supergiant in the SMC. These stars suffer dramatic physical 

' changes on timescales of months - when they are at their warmest they are also 

. brighter, more luminous, and show an increased amount of extinction. This 

L| ■ variable extinction is characteristic of the effects of circumstellar dust, and can 

.£h ! be connected with sporadic dust production from these stars in their coolest 

. states. We suggest that these unusual properties are indicative of an unstable 

' (and short-lived) evolutionary phase not previously associated with RSGs, and 

consider the implications such behavior could have for our understanding of the 
latest stages of massive star evolution in low-metallicity environments. 



1. Introduction 

Red supergiants (RSGs) are a He-burning phase in the evolution of moderately 
high mass stars (1O-25M0). Until recently, the location of RSGs on the H-R 
diagram was not accurately reproduced by stellar evolutionary theory, with the 
stars characterized as too cool and too luminous to agree with the predictions of 
the evolutionary tracks. This places them at odds with the restrictions imposed 
by the Hayashi limit, which denotes the largest radius a star of a given mass ca n 
have while still remaining in hydrostatic equilibrium ( Havashi HoshU Il96ll ). 
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Stars which are cooler than the tracks allow inhabit the region to the right 
of this limit on the H-R diagram, suggesting that they are not in hydrostatic 
equilibrium. 

There are many explanations as to why RSGs might occupy this "forbidden" 
region of the H-R diagram. RSGs present significant challenges to evolutionary 
theory. The velocities of the convective layers are ne arly sonic and even super- 
sonic in the atmospheric layers, giving rise to shocks ( Frevtag. Steffen. Dorchl 
2002) and invalidating mixing- length assumptions. There is poor knowledge 
of RSG molecular opacities. Finally, the highly-extended atmospheres of these 
stars differ from the plane-parallel geometry assumption of the evolutionary 
models. As a result of these uncertainties, the error could lie in the position 
of the evolutionary tracks on the H-R diagram. Another possibility is that the 
"observed" location of RSGs in the H-R diagram is incorrect as a result of inaccu- 
rate effective temperatures - past effective temperature scales for red supergiants 
have been based on temperatures from broadband colors, which ar e highly sensi- 
tive to surface gravit y and the adopted reddening for the sample ( Massevlll998l ; 
Levesque et al .1120051 ). Finally, it is possible that some stars might truly occupy 



this region, suggesting that they are unusual stars in a unique and unstable 
evolutionary state. 



2. Milky Way Red Supergiants 



To shed some light on this debate, Levesque et al.l ( 20051 ) redetermined the effec- 
tive temperature scale for RSGs in the Milky Way. We fit moderate-resolution 
optical spectrophotometry of 74 Ga lactic RSGs with the new generation o f 



MARCS stellar atmosphere models ( Gustafsson et al. 20031 ; Plez et al 



1992) 



which include a much improved treatment of molecular opacities ( Plez 20031 ; 
Gustafsson et al.ll2003l ). With these models we were able to use the rich and 
temperature-sensitive TiO molecular bands to measure their effective tempera- 
tures. The newly derived physical parameters brought the location of the RSGs 
into much better agreement with the predictions of stellar evolutionary theory 
- there is now excellent agreement between Milky Way supergiants and the evo- 
lutionary tracks. This suggests that, on the whole, the theoretical models do 
a good job of accurately reproducing this stage of massive star evolution in 
the Milky Way. However, there are occasional exception to this that must be 
considered. 



2.1. VY Canis Majoris 

VY Canis Majoris is a peculiar late-type M supergiant in the Milky Way. It has 
several physical characteristics unique among its counterparts, mostly related to 
its unusual circumstellar environment. It has a large IR excess, making it one 
of the brightest 5-20 fim objects in the sky, indicative of a dust shell or disk 
being heated by the star ( Herbig 19701 ). Such a dust shell is overwhelmingly 
present - the RSG is surrounde d by an asymmetric dust-reflection nebula that 



extends 8"-10" from the sta r (jMonnier et al.1 Il999l ; ISmith et all l200ll : ISmithl 



20041 ; iHumphrevs etal.1 1200511. Its mass -loss rate is quite high for a RSG, about 



2 x 10 4 M© yr 1 ( Danchi et al.l [i~994l l . The photometric history of this star 
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dates back to 1801; the star has faded by about 2 mag during that time bu t has 

showed little change in its color ( Robinsonlll970l . [l971 ; iMassev et al.ll2006l). 

This star is often assumed to have extraor d inary physical propertie s ( Le Sidaner &: Le Bertrei 
19961 : ISmith et al.1 l200ll : iMonnier et aD 120041 : iHumphrevs et al.l I2Q05I V Its as- 



sumed effective t emperature for many ye a rs wa s a remarkably cold 2800 K, a 
value adopted by iLe Sidaner &: Le Bertrei ( 19961 ) said to be based on the star' s 
assigned M4 I spectral type (the effective temperature scale of lDvck et al.l ( 19741 ) 
is cited, but the 2800 K effective tempearature they report is much cooler than 
the 3400 K that the Dyck scale associates with M4 supergiants). VY CMa's 
luminosity is often cited as approximately M\, \ = -8.5 to -9.5. With such a 
high luminosity and cold temperature, this places VY CMa in a considerably 
cooler and more luminous position that current evolutionary tracks allow, in- 
habiting the previously mentioned Hayashi forbidden region (Figure 2).. Such 
a star would have a radi us of 1800 to 30 R^ , making it far bigger than the 
largest RSGs described in lLevesque et all ( 20051 ) and far outside the limitations 
of evolutionary theory. 

A proper understanding of this star's physical pr operties hinges on an ac- 



curate determination of its effective temperature. In Massev et al. ( 20061 ) we 



used the same method of model fitting to redetermine VY CMa's temperature 
spectroscopically. A satisfactory fit of T e fj = 3650 K was quickly obtained. We 
reclassified the star's spectral type as an M2.5 I based on the strength of the 
TiO bands typical l y use d for such a classification (AA 6158, 6658, 7054; see 
iJaschek k, Jaschekl 19901 ) - this spectral type and effective temperature agree 
excellently with the effective temperature scale for Milky Way RS Gs. Since 
many earlier studies describe VY C Ma as an M3 I to M5 I superg iant (jJovill942l : 
IWallersteinlll958l ; iHumphr evsll 1 9 74l ) based on blue (photographic) spectrograms, 
we also used the TiO bands further to the blue (AA 4761, 4954, 5167) to assign 
a spectral type and temperature, and arrived at a spectral type of M4 I (in 
good agreement with earlier studies) and a T e g- = 3450 K, a cooler (and poorer 
overall) temperat ure that is still much hig h er tha n the commonly-adopted 2800 
K temperature of ILe Sidaner Sz Le Bertrei ( 19961 ). From these temperatures we 
calculated a M^ Q \ = -7.2 f or the red TiO band s, and M^ Q \ = -6.9 for the blue 
TiO bands, as described in lMassev et al. ( 20061 ). 

With these new physical parameters, we found that VY CMa now occupied 
a much less surprising region of the H-R diagram (see Figure 1). The new T e g 
and Mboi derived from spectral fitting place the star comfortably in agreement 
with the rightmost reaches of the 15 M evolutionary tracks. 

However, our method seems to have underestimated the luminosity of the 
star, and Roberta Humphreys (private communication) was kind enough to call 
this to our atten tion. We can d emon strate the problem like so: 

From JHK, IMassev et aD (12005T ) were able to determine the properties of 
the dust shell, and found an effective temperature of 760 K and an area that is 
2155 times that of the star's area. This leads to an interesting contradiction: the 
luminosity of the dust shell must be 4x that of the star, or about 2.3xlO 5 L0. If 
the star's luminosity is the only source of the heating of the dust, and the dust 
is in thermal/radiative equilibrium with its surrounds, then this is very hard to 
understand. 
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Figure 1 . Placement of VY CMa in the H- R diagram. The fille d circle cor- 
responds to the parameters of VY CMa from iSmith et al.l (|2001l ) . The stars 
represent the results of our work, based on fitting the MARCS models to our 
spectrophotometry - the five-pointed star corresponds to the parameters de- 
rived from fitting the red TiO bands, while the four-pointed star corresponds 
to the paramet ers derived from fitt i ng th e blue TiO bands. The evolutionary 
tracks are from lMevnet fc Maederl (j2003) and correspond to models with no 
initial rotation (solid lines) and 300 km s _1 of initial rotation (dashed lines). 
The solid diagonal lines at the upper right represent lines of constant radii. 
The Hayashi limit is represented by the dotted vertical line. 

It is clear that something has caused us to underestimate the luminosity 
of the star. One explanation would be if the star's dust had a larger grain 
size distribution than normal, and that its extinction was much "greyer" than 
normal as a result. This would require a lot of extra extinction - about 1.5-2.0 
mag - but would present one way out of this conundrum. 

Despite our error in determining the luminosity of the central star in this 
unusual object, it should be pointed out that a change in bolometric luminosity 
on the H-R digram shifts the location of VY CMa upward on the H-R diagram 



Unusual Magellanic Cloud Red Supergiants 



5 



but does not disrupt its position with respect to the Hayashi limit, as its temper- 
ature still places it solidly at the edge of higher-mass evolutionary tracks, where 
it remains in hydrostatic equilibrium. Even if the proposed upper limit on the 
bolometric luminosity of Mboi = -9.5 is considered, we can still see that the star 
would simply move from being a stable highly-evolved 15 M RSG to being a 
stable highly-evolved 25 M RSG - its maximum radius would still only about 
about 1200 R© and it would still remain comfortably in hydrostatic equilibrium 
at the limits of the Hayashi track. 



3. Magellanic Cloud Red Supergiants 



With this newfound agreement between the Milky Way RSGs and their evolu- 
tionary tracks, we turned out attention to a si milar problem fo r RSG s in the 
Magellanic Clouds ( Levesque et al.ll2006l ). From iMassev Olsenl 1)20031 ) we find 
that again the evolutionary tracks for the Large and Small Magellanic Clouds 
(LMC and SMC, respectively) do not extend to cool enough temperatures to 
accomodate the "observed" locations of RSGs i n the H-R diagram. W e applied 
the same fitting method described above and in Levesque et al. ( 20051 ). obtain- 
ing moderate-resolution spectrophotometry of 36 LMC RSGs and 37 SMC RSGs 
using the R-C spectrograph on the CTIO Blanco 4 m telescope in November and 
December of 2004. We fit these data with the MARCS models to redetermine 
the effective temperature scales for RSGs in the Clouds. The new effective tem- 
perature scales (see Figure 2) brought the MC RSGs into much better agreement 
with stellar evolutionary theory. The current agreement in the LMC is good, 
while the agreement in the SMC is improved but still not satisfactory, with SMC 
RSGs showing a considerably larger spread in effective temperatures across a 
given luminosity than their LMC counterparts. Such a spread is expected, how- 
ever, due to the lar ger effects of rotational mixing in lower metallicity stars - for 
more discussion see Levesque et al.l ( 20061 ). 



4. Late- Type Red Supergiants in the Clouds 



When comparing RSGs in the Milky Way and the Clouds, lElias et all ( 19851 ) 



note an interesting shift in the spectral types of these stars, with the average 
RSG spectral subtype shifting towards earlier types at lower metallicities. More 
precisely, the average RSG subtype is fo und to be K5-K7 I in t h e SMC, Ml I in 
the LM C, and M2 I in the Milky Way (IMassev Olsenl l2003h . iLevesaue etaf 



(2006) find two distinct explanations for this shift. The spectral subtype of late 



K- and M-type stars is primarily based on the strengths of the TiO bands, which 
are sensitive to both effective temperature and chemical abundances. We see 
that the average spectral subtype shifts to earlier types (and weaker TiO) bands 
in t he lower-metallici ty environments (Z/Zq = 0.2 for the SMC and Z/Z & = 
0.5; IWesterlundl [ l9 97), suggesting that this shift in type is symptomatic of lower 



chemical abundances and weaker TiO lines at a given effective temperature. This 
is well- illustrated in Figure 2, where we can see that a RSG with an effective 
temperature of 3650 K would be assigned a spectral type of M2 I in the Milky 
Way, Ml. 5 I in the LMC, and K5-M0 I in the SMC. It is important, however, 
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Figure 2. Effective temperature scale for the Milky Way (s olid line), LMC 
(dashe d line), and SMC (dotted-dashed line), adapted from iLevesaue et all 
(2006). The three scales together allow us to probe the effects of metallicity 
on spectral type. For example, a 3650 K star would be an M2 1 star in the 
Milky Way, Ml. 5 1 in the LMC, and K5-M0 I in the SMC. 



to note that this explanation does not accomodate the full range of effects that 
metallicity impo ses on massive stellar evolution. 

Figure 10 of lLevesaue et al.l (|2006l ) compares the evolutionary tracks for the 
Milky Way, the Z/Z & = 0.5 LMC, and the Z/Z & 0.2 SMC. There is a clear shift 
of the Hayashi limit - the coolest tip of the tracks - to warmer temperatures 
at lower metallicity. For a 15-25 M star we see a shift in the coolest effective 
temperatures for RSGs of about 100-150 K from the Milky Way to the LMC, and 
about 500 K from the Milky Way to the SMC. In agreement with the speculation 
offered bv lElias et al.l ( 19851 ). this does an excellent job of explaining the observed 
shift in spectral subtypes. Following this explanation, one would expect that the 
shifting of the Hayashi limit to warmer tempeartures at lower metallicity would 
therefore impose a hard limit on how cold, and hence how late-type, RSGs can 
be in a particular environment. 

Despite this expected restriction imposed by the Hayashi limit, there are 
several late-type RSGs in the Clouds that are considerably later than the av- 
erage spectral subtype, cooler than the current evolutionary tracks allow and 
occupying the "forbidden" region to the right of the Hayashi track. We spec- 
ulated that these might not all be true members of the Clouds - instead, they 
could be foreground halo giants. Alternatively, they could present a challenge 
to our understanding of evolutionary theory. We decided to investigate these 
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most unusual Cloud RSGs in more detail ( Levesque et al.l [20071 ) . We obtained 



moderate-resolution spectrophotometry of seven late-type RSGs in the LMC 
and five late-type RSGs in the SMC in December of 2005, using the RC spec- 
trograph on the 4-m Blanco telescope. We fit the spectra using the MARCS 
stellar atmosphere models described above and probed their physical properties 
to better understand the nature of these unusually cool stars. The stars selected 
in our sample, and the physical properties derived from the spectral fitting, 
are gi ven in Table 1 - s evera l stars in the sample had been previously observed 
in the iLevesque et ah ( 20061 ) November/December 2004 run but had not been 
published. 



4.1. Membership in the Magellanic Clouds 

One would expect that a magnitude- and color-selected sample of RSG candi- 
dates could potentially be contaminated by foreground dwarfs or halo dwarfs 
and giants. Of these three, foreground dwar fs are by far t he major contaminant 



by 

in studying RSGs in Local Group galaxies (|Massevl 119981). and can fortunate! 



Massev & Olsenl (200 



3) 



be easily recognized on the basis of radial velocities 
obtained precision radial velocities for a sample of red stars seen towards the 
Clouds using the CTIO 4-m telescope, and found that most of these stars had ra- 
dial velocities consistent with those of the Magellanic Clouds, excluding a small 
fraction (11% for the SMC and 5.3% for the LMC) which had much smaller 
radial velocities and could be readily identified as foreground dwarfs. The re- 
maining stars could be tentatively identified as true RSG members of the Clouds 
- however, most of the apparent radial velocity of the Clouds is simply a reflec- 
tion of the sun's motion, so t his RSG sample could be contaminated by red 
stars in the Milky Way's halo. ILevesque et al.l ( 2006') estimate that this would 
be a few percent or less, but we reconsidered the issue when examining these 
late-type Cloud RSGs. 

The sample of stars in ILevesque et all ( 20061 ) and ILevesque et al.1 ( 20071 ) 
mostly have 12 < V < 14, with a few fainter stars. Their B — V colors are greater 
than 1 .6 (see table 2). According to an updated version of the Bahcall Soneiral 



(1980) model provided by Heather Morrison, we expect a surface density of halo 
giants of about 0.2 ± 0.15 deg~ 2 in this magnitu d e/ col or range towards either 
the LMC or the SMC. The area of the iMassevI (|2002l ) survey was 14.2 deg 2 
towards the LMC, and 7.2 deg 2 towards the SMC, so we might expect a 0.6% 
contamination by halo towards the LMC, and a 0.9% contamination towards 
the SMC. We therefore expect only a fraction of a st a r in o ur entire Cloud 
RSG cample of 85 st ars - 73 stars from ILevesque et al. ( 20061 ) and 12 stars in 
ILevesque et al.l (|2007l ). For the sample of 12 late-type RSGs, a 1% contamination 
(a tenth of a star) is likely a large overestimate, given that the vast majority 
of halo contaminants have B — V < 1.8 while all but two of our late-type stars 
have B — V > 1.8. 

Finally, a more precise kinematic test can be applied to determine if the 
LMC stars discussed here follow the radial velocities of other RSGs as a function 
of spatial position in the LMC. The kinematics of the SMC are quite complex, 
and so we restric t this argument to the L MC, where the kinematics are relatively 
we ll understood (Olsen Massevll2007h . A histogram of RSG velocities derived 
bv lOlsen fc Massevl (120071 ) for the LMC (see Figure 1 of lLevesaue et alJl2007l ) 
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demonstrates that the late-type LMC RSGs in our sample follow the kinematics 
of the galaxy, behavior that we would not expec t in the case o f halo giant 
contaminants - for a more detailed discussion see iLevesque et al. ( 20071 ). We 



therefore conclude that our sample of late-type stars examined here is unlikely 
to contain foreground objects, and are in fact true Cloud RSGs. 



4.2. Photometric Variability 

We found that all 12 of the objects in our late-typesample demonstrated large 
variability in their V magnitu des. Using photo metry from the All Sky Auto- 
mated Survey (ASA S) project ( Pojmanskill2002l ) along with the CCD photome- 
try of IMassevI (|2002l ) and val u es der ived from our own spectrophot ometry using 
the V band curve of IBessell ( 19901 ) and the zero-points given by iBessell et al.l 
(1998),, we calculated a AV^ for each of our stars in bot h the 2004 and 200 5 
samples. Normal RSGs are known to be variable in V ( Josselin et al.l 12000) . 
and a sample of seventy RSGs from ILevesque et al.l ( 20061 ) show an average Ay 
of 0.9 mag. Each of the late-type stars described here show larger variations, 
averaging a Ay of 1.3 for the SMC sample and 1.6 for the LMC sample. The 
most extreme case of V variability is found in the case of SMC 050028, more 
commonly known as HV 11423, which varies by nearly 2 mag. Many of the 
late-type RSG lightcurves also seem to suggest some level of quasi-periodicity. 

In cont rast, we find much lower variability at K, which is also typical of 
other RSGs (I Josselin et all 200Cft: comparing the values from 2MASS and DENIS 
( Kimeswenger et al J 12004 ). we find an average difference in K between the two 
surveys of about 0.08 mag. 

Finally, we find that five of our stars are unusually variable in B — V, with 
the 2005 spectrophoto metry di f fering by several tenths of a magnitude from 
the 1999/2000 values of IMassevI (12002) B - V is not very sensitive to effective 
temperature in RSGs (see Massev et al.ll2007l ). and we suggest that the variations 
we see here are instead indicative of changes in the am ount of circumstella r 
dust, causing differences in the reddening. As described bv lMassev et al.l ((2007), 
dust around RSGs results in a significant amount of circumstellar extinction, 
amounting to several magnitudes in extreme cases. These variations suggest 
ep isodic dust ejection on the timescale of a few years, consistent with the study 



of iDanchi etaD (|l994l ) 



It is important to note that these dust-related changes in B — V do not 
correlate simply with changes in the V magnitudes. This shows that the V- 
band variability is not merely symptomatic of episodic dust ejection, but is in 
fact a result of true physical changes in the star. The time-resolved observations 
of these stars, including their V magnitudes, are given in Table 2. 



4.3. Spectral Variability 

We were intrigued by the large discrepancies in spectral subtypes assigned to 
some of our stars by previous studies and our own past work. While determin- 
ing spectral types involves a small degree of subj ectivity, there was no similar 

/p es and 

we determined for the 2004 observations from ILevesque et all (2006). 



disagreement between the iMassev Sz Olsenl (120031) spectral types and the ones 

Indeed, 



spectral variability of a type or more is unheard of in RSGs. 
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This disparity in spectral types is most profoundly illustrated in the case 
of HV 11423, which also stands out as the RSG that is the most highly variable 
in V. When we observed the star in 2004 we assigned it an unusually early 
spectral type of K0-1 I, corresponding to a temperature of 4300 K. When it was 
reobserved in 2005 its spectrum had changed dramatically, yielding a spectral 
type of M4 I and a temperature of 3500 K. An additional spectrum of HV 11423 
was obtained in September of 2006 with the CTIO 1.5-m telescope, thanks to the 
courtesy of F. Walters and the SMARTS observing queue. While this spectrum 
did not have sufficient wavelength coverage to determine a spectral type or 
effective temperature, it agreed excellently with the K0-1 I spectrum from 2004. 
Furthermore, a December 2001 archival spectrum of HV 11423 from the Very 
Large Telescope corresponds to an even later (and hence, cooler) spectral type 
of M4.5-5 I, and HV 11423 was assigned a spectral type in M0 I in both October 
1978 and October 1979. The M4-5 I spectral types are by far the latest type seen 
for an SMC supergiant, and the corresponding temperatures in these late- type 
states place HV 11423 well outside the limits of hydrostatic equilibrium. 

We also find three other stars in our sample that exhibit evidence of spec- 
tral variability between 2004 and 2005: SMC 046662 (M2 I to K2-3 I), SMC 
055188 (M2 I to M4.5 I ) LM C 170452 (M4.5-5 I to Ml. 5 I) - from compar- 
ing the iMassev fc Olsenl (|2003l ) spectra observed in 2001 to our 2004 and 2005 
observations, we can confirm that these changes are truly symptomatic of differ- 
ences in line strength, and therefore changes in the stars' physical properties (for 
comparison, LMC 148035 was found to vary between M4 I and M2.5 I between 
2001 and 2005, but in this case there was no compelling difference in the spectra 
and we concluded that the disagreement was a result of lMassev &: Olsenl (|2003l ) 
assigning too late a spectral type). All four of these stars - HV 11423, SMC 
046662, SMC 055188, and LMC 170452 - undergo substantial changes in their 
effective temperature over the course of a single year, a phenomenon that has 
not previously been observed in RSGs. For the three SMC stars, we can also 
conclude that when the stars are at their hottest they are also more luminous 
(the 2004 spectrum of LMC 170452 has incomplete spectral coverage, ranging 
from 4000A to 6500A as compared to the 4000A - 9000A coverage of the rest 
of our stars; the data in this case is sufficient for assigning a spectral type but 
insufficient for calculating physical properties such as T e g- and Ay). Finally, we 
are able to estimate the amount of extinction for the three SMC stars, and find 
the Ay is larger when the stars are hotter. This reddening is higher than what 
is typically seen for OB stars in the Clouds, further supporting the argument 
that these effects are symptomatic of sporadic circumstellar dust production. 
Both HV 11423 and S MC 055188 are am ong the only four known SMC RSGs 
that are IRAS sources (|Massev et al.ll2007T ). suggesting that we are m fact seeing 
thermal emission from these stars' circumstellar dust. These variations in the 
stars' physical properties are given in Table 2. Changes in the spectra of the 
unusually variable SMC stars between 2004 and 2005 are illustrated in Figure 
3. 



4.4. Placement on the H-R Diagram 



In Figure 4 we place our full sample of LMC and SMC late- type RSGs on the H-R 
diagram, accompanied by stellar evolutionary tracks of appropriate metallicity. 
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Figure 3. Comparisons of our three variable SMC RSGs' spectra in 2004 
(left) and (2005) right, for HV 11423 (top), SMC 046662 (middle), and SMC 
055188 (bottom). The data are shown as a solid black line, with the fits by 
the reddened MARCS models overlaid as a dotted gray line and the model 
parameters are given in the lower left. 



It is clear that stellar evolutionary theory is not in agreement with the observed 
parameters of these late-type stars, with all of the stars lying on or to the 
right of the Hayashi limit. The location of the RSGs as derived from spectral 
fitting is, on average, 275 K coller than the tracks allow for the LMC and 541 
K cooler for the SMC. When Milky Way RSGs such as VY CMa occupied this 
forbidden region of the H-R diagram, the disagreement was rectified by a more 
accurate determination of the stars' effective temperature or by a reevaluation 
of other physical properties. By contrast, it is evident in this case that the 
location of these evolutionary tracks does not accommodate the full range of 
RSG properties and behaviors in these lower-metallicity environments. This 
is further highlighted by the fact that the discrepancy is worse in the case of 
the SMC, where the lower metallicity will enhance the effects of rotation on 
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the luminosity of the evolutionary tracks (jMaeder Mevnetl l200ll ) due to the 
effects of mixing. Even the high-rotation (300 km s _1 ) tracks are not sufficient 
to reproduce the location of these stars, showing that this increased discrepancy 
in the SMC remains unexplained. 
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Figure 4. Location of our sample of late-type MC RSGs on the H-R diagram 
as compared to the evolutionary. We show the LMC (left) and SMC (right) 
RSGs with the physical parameters derived from fitting of the MARCS models 
to the spectrophotometry. We also show evolutionary tracks of the appropri- 
ate metallicity. The non-rotation evolutio nary tracks are s hown as solid black 
lines and come from iSchaerer etHI (|1993[ ) for the LMC and lMaeder fc Mevnetl 
(f200lh for the SMC. The new models, which incl ude 300 km s 1 ini t ial ro - 
tation, are shown as dashed lines and com e from iMevnet fc Maederl (|2005l ) 
for the LMC and iMaeder fc Mevnetl (|200lD for the SMC. The Hayashi limit 
is represented by the vertical dotted line. In the case of the SMC, we use 
symbols to represent the 2004 (open) and 2005 (filled) observations of our 
three variable SMC stars: HV 11423 (triangle), SMC 046662 (square), and 
SMC 055188 (pentagon). 



We also compare H-R diagrams of the late- type RSGs in the SMC in 2004 
and 2005, illustrating the change in position of HV 11423, SMC 046662, and 
SMC 055188. Surprisingly, we see that HV 11423 actually moves out of hydro- 
static equilibrium between 2004, when it is positioned at the rightmost reaches 
of the evolutionary tracks, and 2005, when it has moved well into the forbidden 
region to the right of the Hayashi track. All three of the stars undergo these 
unusual physical variations while in the Hayashi forbidden region. 



5. What's Going On? 



HV 11423, SMC 046662, SMC 055188, and LMC 170452 all exhibit remarkably 
similar behavior: cool stars, inhabiting the forbidden region to the right of the 
Hayashi track, that show large variability in spectral type, V magnitudes, and 
extinction, presumed to be from episodic production of circumstellar dust. These 
stars suffer changes in effective temperature and bolometric luminosity on the 
timescales of months; in every case, when they are at their hottest they are 
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also brighter, dustier, and more luminous. These stars all occupy the Hayashi 
forbidden region of the H-R diagram, where we would not expect stars to be in 
hydrostatic equilibrium. One would expect hydrostatic instability to lead to this 
sort of variability and behavior, but the precise physical explanation for these 
phenomena and the mechanisms which may lead to such an instability remains 
a mystery. 

It is possible that such behavior is simply of extremes, with the 

coldest stars in any environment moving to the right of the Hayashi track and 
undergoing these unique physical variations. Conversely, it could be possible 
that such extremes are only attainable in low-metallicity environments where 
the Hayashi track moves to warmer temperatures. Similarly, such behavior could 
simply "ignite" only in low-metallicity environments, a function of the complex 
interrelation between massive stellar evolution and the interstellar medium envi- 
ronment. The question of how such environmental factors might correlate with 
this behavior can be solved by observations of the latest-type Milky Way red 
supergiants, as well as observations of RSGs in other low-metallicity environ- 
ments and continued monitoring of these unusual stars as a means of obtaining 
a more complete understanding of the physical phenomena being observed. Bet- 
ter knowledge of this unusual behavior could shed new light on the nature of 
late-type massive star evolution, particularly at low metallicities, and highlight 
properties and environmental effects that the next generation of evolutionary 
tracks should strive to accomodate. 
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Tabic 1. Program Stars and Physical Properties 
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"From merged 2004 and 2005 data. 
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Table 2. 
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"Measurements from 2,451,186-2,451,999 were obtained using the CCD images described in 
iMassevI (|2002f) . with uncertainties of 0.01 mag or smaller. Later measurements come from our 
2004 and 2005 spectrophotometry. 

6 The V magnitudes quoted for the 2,452,188 HJD's comes from ASAS observations on the 
closest approximate date to the observations, from 2,452,184-2,452,187. 
c From merged 2004 and 2005 data. 



